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SECTION  I 
INTRODUCTION 

The  active  cooling  of  re-entry  nose  tips  may  be  the  only  practical 
method  of  maintaining  a shape  stable  nose  in  a severe  weather  environment. 
The  coolant  provides  protection  from  the  thermal  and  particle  environ- 
ments, and  the  degree  of  protection  depends  on  how  effectively  the 
available  coolant  is  used.  The  transpiration-cooled  nose  tip  (TCNT)  is 
one  means  of  actively  cooling  the  nose  tip.  By  minimizing  the  coolant 
necessary  for  thermal  protection,  more  coolant  is  available  for  particle 
protection.  One  of  the  options  available  for  reducing  the  amount  of 
coolant  for  thermal  protection  is  to  optimize  the  external  shape  of  the 
nose  tip.  The  objective  of  this  investigation  is  to  determine  the 
external  shape  of  a TCNT  that  requires  the  minimum  amount  of  coolant  over 
a re-entry  trajectory. 

The  general  problem  of  minimizing  the  aerodynamic  heating  to  a 
hypersonic  vehicle  has  been  considered  by  various  authors  (References  1, 

2,  3,  and  4).  All  of  these  studies  used  the  calculus  of  variations 
method  and  assumed  a modified  Newtonian  pressure  distribution,  either 
completely  laminar  or  completely  turbulent  flow,  and  constant  freestream 
conditions.  The  specific  case  of  the  TCNT  was  considered  by  Baker  and 
Kramer  (Reference  5).  They  also  assumed  a modified  Newtonian  pressure 
distribution,  and  either  completely  laminar  or  completely  turbulent  flow. 
Instead  of  constant  freestream  conditions,  however,  their  trajectory 
parameter  beta,  WT/CpAg,  was  assumed  to  be  a constant.  All  but  one  of 
these  studies  concluded  that  the  flat  face  - sharp  corner  shape  minimized 
the  total  heat  load  to  the  vehicle,  and  that  for  a given  fineness  ratio 
there  was  an  optimum  flat  face  height.  Furey's  (Reference  3)  results 
differed  in  that  the  flat  face  had  a small  "bubble"  at  the  stagnation 
point.  Hull  (Reference  6)  has  pointed  out  the  mathematical  restrictions 
that  result  from  using  the  above  assumptions  with  the  calculus  of 
variations  method.  These  restrictions  make  the  results  in  the  literature 
mathematically  questionable.  Baker  and  Kramer  did  take  Hull's  results  into 
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consideration  in  their  investigation  but  they  assumed  that  the  flat  face 
nose  was  the  optimum  shape  without  ever  proving  It. 

The  results  in  the  literature  may  also  be  questioned  from  a practical 
viewpoint.  The  modified  Newtonian  pressure  distribution  is  of  course 
accurate  for  the  hemisphere  but  not  for  the  blunter  shapes  such  as  the 
flat  face.  The  boundary  layer  flow  over  the  vehicle  is  usually  a 
combination  of  both  laminar  and  turbulent  and,  as  turbulent  flow  is 
associated  with  much  higher  heating,  the  effect  of  the  nose  shape  on 
transition  is  an  important  factor.  The  freestream  conditions,  as  well 
as  beta,  vary  throughout  the  trajectory.  Assuming  that  either  of  these 
is  a constant  removes  the  effect  of  the  nose  shape  on  the  trajectory  time 
which  is  another  factor  effecting  the  total  coolant  required.  In  view  of 
the  above,  it  was  decided  to  reconsider  the  optimum  TCNT  problem  and  use 
the  variation  of  parameters  method  rather  than  the  calculus  of  variations 
method.  The  variation  of  parameters  method  requires  the  calculation  of 
the  total  coolant  needed  by  a given  nose  shape  during  reentry.  The 
optimum  shape  is  the  shape  requiring  the  least  amount  of  coolant.  This 
method  allows  the  optimum  shape  to  be  more  realistically  determined  as 
the  assumptions  mentioned  above  do  not  have  to  be  made. 

As  the  optimum  shape  is  determined  by  comparing  the  amount  of 
coolant  required  by  one  shape  with  that  required  by  another  shape,  only 
relative  accuracy  is  required  in  the  calculations.  The  simpler  empirical 
methods  available  in  the  literature  were  used  to  calculate  the  aerodynamic 
heating  and  drag.  Some  of  these  methods  are  valid  only  for  hypersonic 
speeds,  hence  this  Investigation  is  restricted  to  impact  Mach  numbers 
greater  than  four.  Three  different  families  of  nose  shapes  were  con- 
sidered - the  oblate  ellipsoid,  the  flat  face  - round  shoulder,  and  the 
spherical  arc  - round  shoulder.  These  families  are  bounded  by  the  flat 
face  - sharp  corner  at  one  extreme  and  the  hemisphere  at  the  other 
extreme.  A ten-degree  half-angle  cone  was  used  as  the  aft  body  for  all 
three  shapes.  The  transpirant,  water,  was  assumed  to  cool  only  the  nose 
as  the  aft  body  heating  is  usually  handled  by  an  ablative  material. 

The  wall  temperature  of  the  nose  was  assumed  to  be  2000°R.  At  each 
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altitude  the  coolant  flow  and  total  drag  on  the  vehicle  were  calculated. 
The  trajectory  equations  were  used  to  calculate  the  next  altitude, 
velocity,  reentry  angle  and  range  for  a given  time  step.  A time  step  was 
chosen  such  that  the  drag  coefficient  was  nearly  constant.  The  coolant 
used  over  the  entire  trajectory  was  determined  by  numerically  integrating 
the  amount  of  coolant  used  at  each  time  step. 


The  geometric  considerations  are  presented  In  Section  II.  The 
Invlscld  analysis  is  covered  In  Section  III.  This  Includes  the  shock 
stand  - off  distance  as  well  as  the  properties  behind  the  shock  and  at 
the  body  stagnation  point.  The  viscous  calculations  are  discussed  In 
Section  IV  and  the  trajectory  in  Section  V.  The  validation  of  the 
computer  results  Is  presented  in  Section  VI.  In  this  section,  the 
computer  program  results  are  compared  to  data  in  the  literature  or  values 
obtained  from  hand  calculation.  In  Section  VII  the  results  of  this 
investigation  are  discussed,  and  in  Section  VIII,  the  conclusions  are 
presented. 
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SECTION  II 


VEHICLE  GEOMETRY 


The  equations  used  for  calculating  the  various  geometries  are 
summarized  in  this  section.  The  independent  variable  for  the  nose  is  y 
while  the  arc  length,  s,  Is  used  as  the  independent  variable  on  the  aft 
body. 


1.  OBLATE  ELLIPSOID 


Sketch  1.  Oblate  Ellipsoid 


The  equation  for  an  ellipse  in  the  coordinate  system  shown  above  is 


The  intersection  of  the  nose  with  the  cone  aft  body  occurs  when  9 
The  coordinates  of  this  junction  are 
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The  shoulder  point  location,  ISH,  is  given  by 


The  arc  length  along  an  ellipse  can  be  conveniently  expressed  in  terms 
of  the  elliptic  coordinates  (C »n ) • 


Sketch  2.  Elliptic  Coordinates 


cos  u 

This  is  an  elliptic  integral  of  the  second  kind 


MMl 
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Using  y as  the  Independent  variable,  values  of  x,  0,  n»  are  given  by 


COS  U 


(9a) 

(9b) 

(9c) 


2.  FLAT  FACE  - ROUND  SHOULDER 
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The  Intersection  point  between  the  shoulder  and  flat  face  Is 


The  Intersection  of  the  corner  and  aft  body  Is 


The  shoulder  point  (ISH)  Is 


The  center  of  the  shoulder  (x_.  y.)  is 


The  equations  for  x»  9»  and  s on  the  flat  face  are 


t ;■  * m 

I 


mi  I 


I 


: J 

1 i 
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Intersection  of  arcs  (xj » yj) 


.-1  r *0 


>1  = Sin  ' { F-^—  } 

1 RN  rc 


= rn  TV^l 


X]  ■ Rn  0 - cos  ♦]) 


SI  ’ RN  *1 


Intersection  of  corner  with  flank  (xj,  yj) 


♦j  " t - 9c 


y,  = yrt  + r cos  0 
o c c 


Xj  = X0  - rc  sin  0C 

Sj  = si  + rc  (<t>j  - ♦j ) 


Shoulder  location  (Xjsh*  ylSH^ 


XISH  = xo 


yISH  ■ yI  + <x„  - Xj)  8C 

. , *yISH  ' yj) 

sISH  ■ SJ  + sin  0. 


07a) 


(17b) 


(17c) 


(17d) 


(18a) 


(18b) 


08c) 


(18d) 


(19a) 


(19b) 


(19c) 


The  equations  for  the  dependent  variables  over  the  different  regions  are: 
On  spherical  arc  (y  < y j ) 


♦ a sin"1  (y/RN) 

x " rn  - a/rn  ■ yz 


(20a) 

(20b) 
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9 * | - + 
s = <J> 

On  corner  (y.  < y < y,) 

♦ 3 sin'1  (^-^) 
c 

X 3 X0  - rc  cos  <f> 

9 * \ - ♦ 


i 

(20c) 

(20d) 

(21a) 

(21b) 

(21c) 


s = si  + rc  (♦  " ^ 


On  shoulder  (yj  < y 1 yISH) 

X 

s 


(y  - yj) 

tan  9C 


(y  - -yj) 

sin  ec 


9 


9 


c 


(21dl 


(22a) 

(22b) 

(22c) 


It  Is  also  noteworthy  to  consider  the  case  of  the  sharp  corner.  For 
0j  greater  than  10°,  a discontinuity  will  exist  at  the  corner. 


yJ  = yISH  = a 

XJ  = XISH  = RN  ' 7^  " ^ 


For  0,  less  than  10°,  the  nose  will  flare  smoothly  into  the  aft  body 


* 
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SECTION  III 
INVISCID  CALCULATIONS 

In  this  section  the  methods  used  to  calculate  the  invlscid  local 
flow  properties  are  presented.  Empirical  methods  available  in  the  liter- 
ature were  used.  Real  gas  effects  were  accounted  for  by  determining  an 
effective  specific  heat  ratio. 

1.  FREESTREAM  CONDITIONS 

The  frees tream  conditions  at  a given  altitude  were  determined  from 
equations  given  by  Dommasch,  Sherby  and  Connolly  (Reference  7).  The  sea 
level  values  are  as  follows: 

P$L  = 2116.217  lbf/ft2 

p$L  = 0.002377  slugs/ft2 

T$l  = 518.688  °R 

a^  a 1116.444  ft/sec 

g$L  = 32.174  ft/sec2 

R = 20,855,531.0  ft 
er 

The  distance  above  the  earth's  surface  is  the  geometric  altitude, 
hg.  The  distance  from  the  earth's  center  is  the  absolute  altitude,  ha. 
The  geopotential  altitude,  hg  is  given  by 

h = hr;  hr  < 60,000  ft  (28a) 

g b b — 

Rov>h 

h = ; hr  > 60,000  ft  (28b) 

g ha  G 

In  the  geopotential  layers  where  g Is  essentially  constant,  the  pressure, 
density  and  temperature  are  given  as  follows: 

Isothermal  region  (T  = constant): 


6 = <$i  eG 

(29a) 

a = o1  ee 

(29b) 

13 
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I 


Gradient  Region  (dT/dhg  Is  constant): 


T ' T1  + ad  (hg  - hb>  " T1  + adhg 


« ■ (T/T'j ) ’1/adr 
a - o1  (T/T1 ) -[1/adR  + 1] 


(30a) 

(30b) 

(30c) 


where 


hK  - hn 

b 


dT 


RT 


’d  dh. 


and  hb>  6-j  > o-j , T-|  and  y are  given  In  Table  1 


l/adR 


The  freestream  stagnation  properties  are  calculated  from  the 
relationships  using  the  static  values  of  P^,  T^,  and  and  the  free- 
stream values  of  M and  y. 


VT-[’ 

] . (°R) 

Ola) 

p0„  ■ p.  [ 1 

l^.af 

J ft£ 

(31b) 

\ ■ p“  [ 1 

] 

J ft* 

(31c) 

enthalpy  was  determined 

as  follows 

h = Cn  T 

°co  P “ 

V2 

00 

+ 2Jgc 

(32) 

2.  SHOCK  DENSITY  RATIO 


If  hQ  Is  less  than  864.0  BTU/lbm,  the  isentroplc  gas  equations  can 


be  used  to  determine  the  shock  density  ratio.  For  hQ  greater  than  864.0 

00 

BTU/lbm,  real-gas  effects  must  be  considered  and  an  iteration  procedure 
Is  required. 


14 
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Ideal  Gas  (hQ  < 864  BTU/lbm) 


P2  - Y + 1 

p»  (y-D+/ 


To  [2*£  - (y-D3C(y-1)h£  + 2] 


(Y+l )2  M* 


— a 1 + (m2  _ 1 ) 

P Y+l  ^ 

OO  * 

(Y-l)  m!  + 2 


2 YM:  - (y-1) 


(33a) 


(33b) 


(33c) 


(33d) 


Real  Gas  (hQ  > 864  BTU/lbm) 


The  computational  procedure  for  the  real -gas  case  Is  as  follows: 
(1)  Guess  p2  (Isentroplc  relationships) 


<2>  u2  = p-  u»/p2 


2 (3)  P2  ■ P.  ♦ P.  - P2  4 


< 

£ (4)  h2  = h^  + 1 (u*  - uf) 

-(5)  p2  = p2  (h2,  P2);  Bade  (Reference  19) 


The  Iteration  Is  continued  until  the  percent  change  In  the  density  Is 
less  than  1 x 10~  . 

Knowing  the  density,  pressure  and  enthalpy  behind  the  shock,  a 
temperature  Is  calculated  from  the  perfect  gas  equations.  If  the 
temperature  Is  less  than  2500°R,  the  following  equations  were  used. 
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. i 


0.4 


5500) 


(34a) 


1 + 0.4  (^) 


(5500/T) 


[e5500/T.1]2 


u2 


(34b) 


If  is  greater  than  2500°R,  then  the  temperature  behind  the  shock  is 
calculated  from  Hansen  (Reference  8).  Because  of  its  complexity,  this 
method  is  not  presented  here. 


The  stagnation  conditions  behind  the  shock  are  calculated  as  follows: 


"2  ' V/ 


yRT, 


o 

p 

2 


?o„  = P2  + I p2 


■sM] 


P°2  * P°2  ^H°’  P°2^  BadC  ^ReferenCe  19^ 


(34c) 

(34d) 


T = T (P.  » h_  , p.  );  Hansen  (Reference  8)  (34e) 

0^  o 2 O2  °2 


^eff  = £n  (P2/Po2^/£n  (p2/po^ 


(34f) 

M2  = u2</  Yeff  R T2  (34g) 

The  iteration  is  continued  until  the  change  in  M2  is  less  than  0.001. 


SHOCK  WAVE  CALCULATIONS 


The  shock  stand-off  distance,  radius  of  curvature  and  profile  were 
calculated  from  Krasnov  (Reference  9),  and  James  and  Terry  (Reference  10). 
The  shock  stand-off  distance  and  radius  of  curvature  for  the  three 
different  nose  shapes  are  as  follows. 
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a.  Sphere 


Sketch  7.  Hemisphere  Shock  Stand-off 


Shock  stand-off 


S.  = 0.78125  p 


Radius  of  curvature 


R„  = 0.78125 


b.  FI  at- Face  Sharp  Corner 


f [1.0+0.78125^] 


c.  Flat  Face-Round  Corner 


Sketch  8.  Flat  Face  - Round  Shoulder  Shock  Stand-off 


p 1 

Pi 
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So  ■ f * Sf  * tsf  • V -f 


K *, 

sf  sh 


R. 


(37a) 

(37b) 


d.  Oblate  Ellipsoid 

The  shock  stand-off  distance  for  the  oblate  ellipsoid  was 
adjusted  to  fit  the  trends  given  by  Felderman  et  al  (Reference  11). 

S = s e £ei(b/a)  + e2(b/a)2] 
e f 


(38) 


The  constants  e-j , and  e2  are  calculated  from  and  SQ  5.  A simultaneous 
solution  gives 

e1  = 4tn  {0.025  + 1.1814  p } - in  (Sh) 
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Assume  that  SQ  varies  linearly  with  (b/a),  then 


Se  ■ Sf  - (sf  - Ih)(b/») 

(41a) 

% - V5,  * v^f  - (Sf  - shxb/»>J 

(41b) 

where  Sg  is  given  by  Equation  38. 

e.  Spherical  Arc  - Sharp  Corner 

The  shock  stand-off  distance  and  shock  radius  of  curvature  for 
the  spherical  arc  nose  shapes  are  given  in  this  section.  Results 
presented  in  Krasnov  (Reference  9)  show  that  the  shock  stand-off  varies 
linearly  with  (A/a). 


where 

a * T [,  0-  Vl  °-(f-)Z]  <«> 

Equation  42  Is  assumed  to  be  valid  for  when  the  sonic  point  is  on  the 
corner,  i.e.  0 < (A/a)  < A*.  The  sonic  point  is  assumed  to  be  at  the 
corner  for  0j  >_  45°.  Therefore  A*  * 0.41421,  for  a nose  body  radius  (a) 
of  one  inch. 

For  (A/a)  larger  than  A*,  the  shock  stand-off  distance  must 
decrease  monotonically  to  the  value  for  the  hemisphere.  This  is  accom- 
plished by  matching  a second  order  polynomial  to  a straight  line  through 
(s/a)h. 

(s0/a)h  = CQ  + C1  A + C2  A2  (44a) 

where 

c0  ■ J1  - “ Kf  - 2 V <44b> 


20 
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and. 


C * - — r $ _ ? S 1 
L1  9 f 1 n-* 


A = A - A, 


(44c) 


(44d) 


The  match  point  of  the  polynomial  and  straight  line  is  assumed  to  be 


therefore 


AMP  s 


-2.0  (Co  - Sh) 


. Co  + C1  Ah  ' Sh 
C2  - — 


*-amp  " 2 V Al 


MP 


(45) 


(46a) 


where 


also. 


Ah  - 1.0  - A, 


"SA 


C1  * 2,0  tCo  + ClAh  " V 

^amp-2V 


(46b) 


(46c) 


and. 


b = Sh  ‘ mSA  Ah 


(46d) 


The  polynomial  expression  is  used  from  A*  < A £ (A*  + AMp).  The  linear 
expression  ST  = b + mA  is  used  from  (A*  + A^p)  to  1.0. 

An  expression  for  the  shock  radius  of  curvature  is  given  by  Krosnov  as 


Rs  Bs 

f 


V ” [\  - «,J 


(47) 


where  B Is  given  as  the  value  of  R at  A = 0.45  and  is  approximately 

5 5 . 

— 0 20 

equal  to  1.8  to  2 times  Rs.  . For  this  investigation  B Is  taken  as  -g- 

»h0  s » 

R$h  , which  is  somewhat  larger,  and  the  above  equation  is  used  over  the 

full  range  of  A. 


21 


AFFDL— TR-78-1 47 


(x  T 0.75  ? 1-0.07 

0.75  2 ) 0.2741  + 1.15  Cn  (m£-1  ) 

CD  (M„-D>  L °N  “ J 

/ 0.75\ 

’(•°-237cC  ) 


0.646  exp(-0.237  Cr 
n = — p- 


cS'5 

°N  “ 


0.034 


(50b) 


(50c) 


where  Cn  is  the  drag  coefficient  for  the  different  nose  shapes. 
UN 


4.  PRESSURE  DISTRIBUTION 

The  inviscid  pressure  distribution  over  the  reentry  body  was  con- 
sidered in  three  parts:  the  nose,  after  body  and  base.  There  are 
many  methods  available  for  calculating  the  pressure  distribution 
over  the  nose  of  a hypersonic  vehicle.  These  methods  however  are 
complicated  and  time  consuming.  For  the  purposes  of  this  investigation, 
a simplified  method  was  desired  that  would  give  relatively  accurate 
results  between  shapes  as  well  as  short  computer  times. 


a.  Oblate  Ellipsoid 

Vinokur's  method  (References  12,  13,  and  14)  is  used  for 
calculating  the  pressure  distribution  over  the  oblate  ellipsoids. 
This  method  gives  good  results  for  the  flat  faced  shape  but  is  some 
what  less  accurate  for  the  hemisphere.  The  method  is  developed  in 
detail  in  Reference  15,  and  is  summarized  here  for  completeness. 


The  nondimensional  pressure  (P/P 


(2.0  - k) 


!) 


) is  given  by 


1 .0-n 


.2 

’B 


r-  + n" 


(51) 
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and  n are  elliptic  coordinates  given  by 


where  k 


The  terms  A & B are  constants  given  by 
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(use  sign  which 


The  constants  C-j 


gives  X < £$) 


<5*  * » 


1 10 


'8 


P2  " - P1  f,  - °-40 


— o CR 

P3  = 2(S  + ])  [C7  - C6  C^ 


i C2  > ^4 » ^5 » C7  > snd  Cg  ape 


C,  = -o-2 — + tan  cs  - d1 


-1  t*' 


C =3j_il 
L2  10 


2 


C4  " 15  - 73 


-1  r 


(q  + 1) 


+ (5C  + 1)  3 tan  £ 


-2  t2 


C6 


h ” 5.s 


10(1  + 5£|) 

-Tc 


7 5(5  5g  + 1 ) 


Cg  = 10  es  (3  tan_15s-d2)+ 


16  + 50  £?  + 30^ 


U^  + D 


25 
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The  constants  d-j  and  are 


+ tan 


(57a) 


15£p  , 

d, o-5— p-5— p + 3 tan  1 CR  (57b) 

2 (552+1)  (1  + e|)(55j  + 1)  B 


b.  Flat  Face  - Round  Shoulder 

The  pressure  distribution  on  the  flat  face  - round  shoulder 
nose  was  determined  by  assuming  a velocity  distribution  from  the 
stagnation  point  to  the  sonic  point.  A modified  Newtonian  type  equation 
was  used  from  the  sonic  point  to  the  Prandtl -Meyer  match  point. 

The  velocity  distribution  was  assumed  to  be  linear  from  the 
stagnation  point  to  SL>  i.e.,  u = (du/ds)Sps,  where 


sL  = 0.70  (1.0  - rc/a)s*  (5£ 

Beyond  sL,  the  velocity  is  assumed  to  vary  parabolically  to  the  sonic 
point  S*.  i.e.. 


“ = “L  + (?“)sp  !S"Sl>  * °3  'S'Sl)< 


where 


u*  - u.  - (du/ds)-p  (s*  - s.  ) 

C-  = ^ (60) 

(s*  - sL)2 

The  corresponding  pressure  is  determined  from  the  isentropic  relationships. 


A modified  Newtonian  type  equation  is  used  aft  of  the  sonic  point 


on  the  rounded  corner. 


O 


- Kd  sin2  (0*  - 0)  - JD(0*-  0) 
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where. 


Jp  = He*  * slope  at  the  sonic  po1nt 


(62a) 


(P*  - pJ  - V*  p* 

2 * r p 

sin  0*  o2 


(62b) 


The  value  for  Jp  was  determined  as  follows: 


where 


dp*  dp*  du*  ds* 

d0  du*  ds*  de 


(64a) 


du*  _ /du\ 
ds*  \ds/ 


+ 2 C3  (s*  - sL) 


(64b) 


dp*/du*  is  determined  from  the  isentropic  relationships  and  M = 1.0, 

...2  . / .v  y1' 


d£*  _ d£*  dM_  _ -iu^L  /y+1_\ 
d“*  = dM2  dU*  ' d,  \ 2 / 


therefore 


Yu*  r 


Vs*-si> 


(ff 


This  method  is  used  from  the  sonic  point  to  the  Prandtl -Meyer  match 
point. 

c.  Spherical  Arc 

The  spherical  arc  pressure  distribution  was  calculated  in  a 
manner  similar  to  that  for  the  flat  face  because  a modified  Newtonian 
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pressure  distribution  fails  to  give  the  proper  stagnation  point  velocity 
gradient  and  sonic  point  location.  The  linear  velocity  region  is  given  by 


This  gives  the  same  results  for  the  flat  face  and  hemisphere  as 
the  equation  for  the  flat  - round  case.  The  term 


is  required  in  order  to  prevent  an  inflection  point  occurring  in  the 
velocity  distribution.  The  equations  for  the  velocity  distribution  in 
the  subsonic  region  are  the  same  as  for  the  flat  face  - round  shoulder 
case.  The  method  for  calculating  the  pressure  distribution  aft  of  the 
sonic  point  is  also  the  same.  This  method  gives  results  which  are 
identical  to  the  flat  face  - round  corner  nose  for  both  the  hemisphere 
and  the  flat  face  nose  shapes. 

The  methods  for  calculating  the  pressure  distributions  over  the 
various  nose  shapes  are  used  up  to  the  point  where  a match  is  obtained 
with  the  Prandtl -Meyer  pressure.  This  occurs  when  the  pressure  and 
pressure  gradients  are  equal  to  those  given  by  the  Prandtl-Meyer  theory. 
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} 


For  the  oblate  ellipsoid  nose 


dp  _ dp  dn 
d9  dri  d0 

where 


d0 

. J 

1 +^B 

dn 

,2.2  f 
+ n 

1 -n2 

2r)k 

<\  - Bv12>  . 

h - n2 

d0 

(2  - 

kl 

(«|  + n2)2  y 

1 + 4 

X {2Bv(l  - n 2)(4  + n2)  + (Ay  - Bvn2)U2  + 1)} 


(70a) 


(70b) 


(70c) 


For  the  flat  face  - round  shoulder  and  spherical  arc  noses,  the 
pressure  gradient  is 


§ - Kp  sir  2 {e*  - 9)  * Jp 


(71) 


The  match  point  location  is  determined  using  an  Iteration  procedure. 
Knowing  the  match  point  angle,  0mp,  and  Mach  number,  Mmp,  the  Prandtl 
Meyer  relationship  is  used  to  determine  the  Mach  number  distribution 
over  the  nose.  The  isentropic  relationship  is  then  used  to  calculate 
the  associated  pressure. 


The  effects  of  mass  addition  on  the  invlscid  pressure  over  the 
nose  was  neglected  (Reference  16). 

d.  Aft  Body  Pressure  Distribution 

The  pressure  distribution  on  the  aft  body  can  be  determined  from 

a method  given  by  Krasnov  (Reference  9)  which  takes  into  account  the  nose 

bluntness  effects  on  the  aft  body  distribution.  This  method  Is  good  for 

for  x > Xish  + 2a.  For  the  region  between  Xj  and  X$n  + 2a,  an  exponential 

relationship  was  used  to  Insure  that  no  Inflection  points  occurred  between 

X,  and  X + 2a. 

J s u 


29 


AFFDL-TR-78-147 

P - PMiN  = A0  exp  [A,  (x'/a)  ♦ ^ (X/a)2]  (72) 

where  x'  is  measured  from  Xj,  and  PMjN  Is  from  Krasnov. 

The  constants  AQ,  A1 , and  A2  are  determined  from  the  following 
boundary  conditions: 

1 v = a x*  = 0-  P & (dn/del  are  known  from  nose 

x “j*  X PJ  ldp/  calculations. 

2*  x'  = x'min  : ^dx>  a 0 


These  boundary  conditions  give 


Ao 

= PJ  ' PMIN 

(73a) 

A1 

(dp/dx)j  a 

Ao 

(73b) 

A 

- _A1  “'hin 

H 2 

2-0x'min  * " HIM 

\ » J 

In  order  to  determine  A-j , 
from  (dp/de)j 

(dp/dx)j  must  be  known.  This 

can  be  determined 

* (i)0  fe), 

(74) 

Therefore,  (d0/dx)j  must  be  known  for  each  of  the  different  types  of 
nose  shapes. 
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Id.  Oblate  Ellipsoid 


Sketch  10.  Oblate  Ellipsoid  Nose  Notation 


■ Xj»  yj.  and  0j  are  given  by  equations  2 and  8b,  respectively 
2d.  Rounded  Corners  (Flat  Face  or  Spherical  Arc) 


Sketch  11.  Round  Corner  Notation 
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/ii\  = =1 

W/j  rc  cos  6J 

where  0j  Is  given  by  Equation  11c  or  Equation  21c 

This  equation  is  also  applicable  to  the  spherical  arc  - sharp  corner 
case.  The  term  rc  is  replaced  by  R^,  i.e., 

fc)j  = rn  C0S  0J 
3d.  Flat  Face  - Sharp  Corner 

y 
s 

Sketch  12.  Flat  Face  - Sharp  Shoulder  Notation 
on  the  flank  x = (S  - Sf)  cos  0c>  and  dy/d0  = 0. 

For  this  case  use  (dp/ds)j 

(I),  - {§\  (f)0 - £ - nb; 

The  constants  Aq  and  A-j  can  now  be  calculated.  The  constant  A^,  as  well 
as  the  pressure  aft  of  x = X$h  + 2a,  requires  Information  from  the 
universal  curve. 

Krasnov's  universal  curve  method  utilizes  x as  the  independent 
variable  and  P as  the  dependent  variable,  where 
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1/  0 . 2 « 

(tan  Q^1 


(79b) 


K.  = M„  tan  0,. 
00  c 


(79c) 


F 


CP  ■ (p  - fJ/\  p.  £ 


<79d) 


where  y,.  is  measured  from 


This  curve  (Reference  9,  Page  606)  was  digitized  as  follows: 


1.475 

1.240 

1.200 

1.200 

1.230 

1.350 

1.480 

1.550 

1.600 

1.650 

1.680 

1.740 

1.770 

1.780 

1.790 

1.800 

0.13,  with  PM 


The  minimum  was  estimated  to  occur  at  X„ 


1.20. 


Using  Equation  79a, 


XkMIN  ' 2>  *MIN  X5 
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Sketch  13.  Shoulder  Region  Notation 


XMIN  = XSH  + xkMIN 
XMIN  = ^XSH  " XiP  + xkMIN 
^MIN  = xMIN/a 

and  A,  3 - A,  a/(2.0  X'mtw) 


(81a) 


(81b) 


(81c) 

( 81  d ) 


The  method  used  In  computer  program  Is  to  use  the  expotentlal  relation- 
ship, Equation  72,  between  Xj  aid  x^j^»  and  then  for  x > > the 

universal  curve  is  used.  Linear  Interpolation  is  used  for  values  of  x 
that  fall  between  the  tabulated  values.  The  pressure  distribution 
obtained  in  this  manner  was  considered  to  be  adequate  for  the  purposes 
of  this  investigation,  i.e.,  no  corrections  were  made  for  curved  shock 
effects  or  boundary  layer  effects  on  the  aft  body. 


S3 'T  rTMr:i  "T~rr" 
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Knowing  the  universal  pressure  distribution  over  the  body,  the 
value  for  velocity,  temperature,  and  density  can  be  obtained  from 
the  Isentropic  relationships: 

[>  + ^ H2]  * ('"oj)  ' Y <82a> 

n2)  ' (82b) 

*02 


°2 

£-■  «2]  ^ 

Oo 


(82c) 


(82d) 


I 

i 


I 

I 


I 


where  y is  obtained  from  Equation  34f,  and  a 


The  enthalpy  is  obtained  from 


h = H - -5— 

o 2 gcJ 


(82e) 


5.  SONIC  POINT  LOCATION  AND  STAGNATION  POINT  VELOCITY  GRADIENT 

In  order  to  calculate  the  pressure  distribution  on  the  flat  face  - 
round  shoulder  noses  and  the  spherical  arc  noses,  the  stagnation  point 
velocity  gradient  and  sonic  point  location  must  be  known. 

a.  Flat  Face  - Round  Shoulder 

The  sonic  point  location  for  the  flat  face  - round  shoulder  nose 
shapes  was  assumed  to  vary  linearly  with  respect  to  (rc/a)  as  follows: 

s*  ' s*flat  ' ' s*hem(]  <rc/a)  (83> 
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; 1 

♦ 


i 


b.  Spherical  Arc  - Sharp  Shoulder 

For  the  spherical  arc  - sharp  corner  nose  shape,  the  sonic  point 


was  assumed  to  be  fixed  at  the  corner  for  9corner  1 45°.  For  9corner  < 45°, 


the  sonic  point  was  assumed 

to  vary  as  follows: 

s*  ' s*hem1  * b 

1 <Z-  W *b2  <Z- 

■ z )Z 

4hemi ' 

(84a) 

where 

Z = (a/RN)_1  ; 

Zhemf  " '’0 

(84b) 

bl  = a sin"’  j 

!VRt! 

(84c) 

- S9=45°  • 

Shemi  ' bl  {Z9=45° 

Zhemi ^ 

(84d) 

b2  " 

(Z  z 

u9=45°  hemi ' 

The  above  expression  for  s*  was  based  on 

a modified 

Newtonian 

type  variation  for  s*,  i.e. 

> 

* 

S Mod 

Newtonian 

ra  1 
o> 

t/> 

3 

1 

1 1 

in 

1 

1 «" 

(85) 

c.  Spherical  Arc  - Round  Shoulder 

For  the  rounded  shoulder  case,  the  sonic  point  was  assumed  to 


vary  with  respect  to  (rc/a), 


s*  = sSC  • {sSC  -shemiXrc/a) 


(86) 


where  sS(.  is  the  value  for  the  sharp  corner. 
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The  stagnation  point  velocity  gradients  were  determined  from  results 
given  by  Boison  and  Curtiss  (Reference  7).  The  stagnation  point  velocity 
gradient  for  the  different  noses  was  assumed  to  vary  with  respect  to 
(X*/y*)  as  follows: 


['M.1IWJ-  (fclj 


- ffl 


Alat"  (^)hemi 


(87a) 


where 


(87b) 


v " • 


[',(*),]  - f,(s)„] 


are  determined  from  Vinokur's  method. 
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SECTION  IV 

VISCOUS  CALCULATIONS 

This  section  describes  the  methods  used  to  obtain  quantitative 
information  concerning  the  viscous  flow  effects.  The  viscous  calculations 
are  more  difficult  than  the  inviscid  calculations.  The  two  calculations 
are  of  course  coupled  and  under  certain  circumstances  this  interaction 
cannot  be  ignored.  In  addition  there  are  the  added  complications  of 
mass  addition,  transition  and  real  gas  effects.  A solution  to  the 
entire  viscous  problem  is  not  required  here  and  use  is  made  of  the  many 
approximations  that  exist  in  the  literature. 

1.  VISCOSITY 

The  viscosity  is  determined  from  an  empirical  correlation  to 
Sutherland's  formula  developed  by  Hansen  and  Heims  (Reference  18). 


y = l ps  T < 6000°R 

(Vs  T - 6000 °R 


(88a) 


where  Sutherland's  law  provides 


= 0.227  x 10  T +T198  6 (fj-sec)  (88b) 


where  T is  in  °R 


, n A 0.023  T 1800  (1,°  " °*125V  ‘6-5 

1,0  1800.0  I1,0  tanh  1.5  + 0.125A..  . 


1.0+  exp 


W- 14-5- 1-5  V 

^ 0.9  + 0.10  A J 


(88c) 


where 


Ay  = log10  (Pa}  ; Pa  is  in  atm 
K$  > 0.04 


(88d) 


J 
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2.  DENSITY 

The  density  is  determined  from  a method  given  by  Bade  (Reference  19). 


where 


range 


X 

h 

P 

P 


P 

1 

1 

1 


31 .9  < 


h 


< 480 


= 0.70  + 0.04  1og1Q  (P/P1 ) 

= 31.9  RTq  = 1080  BTU/lbm 

= 1 atm  = 2116.22  lbf/ft2 
= 0.124  p$L  = 0.00949  lbm/ft3 

; 0.1  atm  < P < 100  atm 


(89a) 

(89b) 

(89c) 

(89d) 

(89e) 


The  above  equation  is  valid  to  within  5%  at  pressures  considerably 
greater  than  100  atm  over  the  enthalpy  range  indicated. 


For  the  range  3.8  £ h/RTQ  < 31. 9,  0.1  < P(atm)  < 100,  XQ  = 0.94 

(In  this  region  p/P  is  nearly  independent  of  pressure.) 


3.  LOCAL  HEAT  TRANSFER  AND  SKIN  FRICTION 

There  are  many  methods  available  for  estimating  the  heat  transfer 
and  skin  friction  on  a flat  plate  in  a uniform  flow.  These  methods  can 
also  be  used  for  a sharp  cone  through  the  use  of  a simple  correction 
factor,  if  it  is  assumed  that  the  local  body  radius  is  much  greater 
than  the  boundary- layer  thickness.  The  essential  characteristic  of  these 
methods  is  a zero  pressure  gradient.  However,  most  shapes  of  interest 
have  a nonzero  pressure  gradient  at  the  boundary  layer  edge,  which  com- 
plicates the  solution.  For  regions  where  the  pressure  gradient  is  not 
too  strong,  the  advantages  of  the  simpler  zero  pressure  gradient  methods 
outweigh  the  increase  in  accuracy  of  the  nonzero  pressure  gradient  methods 


AFFDL-TR-78-147 


m 


r 
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For  this  investigation,  where  only  a relative  accuracy  is  necessary, 
Eckert's  (Reference  20)  "reference  enthalpy"  method  was  used  to  predict 
the  heat  transfer  and  skin  friction  downstream  of  the  stagnation  point. 
This  method  is  not  only  easy  to  apply  but  is  applicable  to  both  laminar 
and  turbulent  flows.  Introducing  properties  evaluated  at  the  reference 
enthalpy  causes  the  effect  of  compressibility  on  the  friction  factor 
to  vanish,  thus  validating  the  use  of  the  incompressible  solutions  for 
both  laminar  and  turbulent  flows. 

As  determined  empirically  by  Eckert 


h*  = 0.5  (hw  + he)  + 0.22  (hr  - he)  (90) 


where  hw  is  the  wall  enthalpy,  hg  is  the  boundary  layer  edge  enthalpy, 
and  hr  is  the  recovery  enthalpy,  given  by 


h 


r 


(91) 


1/2  1/3 

where  the  recovery  factor,  rf,  is  Pr  for  laminar  flow  and  Pr  ' for 
turbulent  flow. 


The  wall  temperature  of  the  nose  was  assumed  to  be  constant  (2000°R), 

thus  h , can  be  determined  from 
w 


where  Cp  was  assumed  a constant. 


CPTw 


(92) 


The  heat  transfer  at  the  stagnation  point  was  calculated  from  the  Fay 
Riddell  equation 


(q)„0  * °-763  ,>r'2/3 

blowing 


(Vh„> 


(93) 


where  Pr  was  assumed  to  have  a value  of  0.71.  A slight  improvement  can 
be  obtained  by  using  expression  given  by  Fox  & Libby  (Reference  21), 

Pr  = 0.525  T0,07 
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For  laminar  heating  downstream  of  the  stagnation  point,  Lees 
(Reference  52)  equations  were  used  in  conjunction  with  the  reference 
enthalpy  method. 


where 


F*(s) 


'"o'  no  blowing  1 ' 

H fp*  1 ruel 

2 L V.J  W y 

’ fs  P*  u*  ue  ? " 

1 — y ds 

II  W u»  J 


(95a) 


I 

L. 


F*(o)  = -J  f-  M 

1 um  D„  U_  \dsj 


(95b) 


The  laminar  skin  friction  was  calculated  from 


Cf  = 0.664  ^zO/Re^ 


where 


C*  = p*u*/PeUe 


and  Res>  is  the  edge  Reynolds  number 

pus 
_ Me  e 


(96a) 


(96b) 


(96c) 


■ 
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The  turbulent  skin  friction  was  calculated  from  an  expression  given 
by  Langanelll  et  al.,  (Reference  23). 


where 


-f-  - 0.0296  Mf  Cj  Re 


1.176 


(97a) 


(97b) 


(97c) 


The  turbulent  Stanton  number  can  be  calculated  using  Reynolds 
Analogy 

Cf  -2/3 

St0  ■ PrT  (9£ 

The  turbulent  Prandtl  number  (PrT)  was  assumed  to  have  a value  of  one. 


The  turhulent  heat  transfer  was  calculated  as  follows: 


(4)  no 


blowing 


peue  (Vh«>  Sto 


The  above  expression  for  heat  transfer  and  skin  friction  were  corrected 
for  diffusion  as  follows: 


‘ [' t(Le6*”^w]  000i,, 


DCIurb  ■ ['  * °-4  Th^y] 


Where  Le  is  the  Lewis  number,  and  was  assumed  to  be  constant 
(Le  = 1.45).  A value  of  0.53  was  assumed  for  8 corresponding  to 
equilibrium  flow. 


(100b) 
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The  enthalpy  of  dissociation  hQ  Is  given  by  Krasnov  (Reference  9). 


T = 273°K 
o 

<j>  =*  0.1;  T > 1000°K 


(101a) 


(101b) 


The  above  equations  were  corrected  for  mass  addition  by  utilizing 
expression  given  in  Timmer,  Arne,  et  al . (Reference  24). 

For  laminar  flow: 


1.0  - 0.68  W°-4B  + 0.08  W°-4B2 


(102a) 


While  for  turbulent  flow: 


(102b) 


The  specific  heat  ratio,  W,  was  assumed  to  have  a constant  value  of 
1.86.  This  is  for  water  vapor  to  air  and  was  given  in  Gold  (Reference  25). 

The  blowing  parameter,  B,  is  defined  as 

B = li,/peuesto  003) 

The  mass  addition  is  assumed  to  be  the  minimum,  i.e., 

m = d/AHc  (104) 

where  AH  is  the  total  enthalpy  rise  of  the  coolant.  For  the  case  of 
c 

water,  a constant  value  of  1150  BTU/lbm  was  assumed. 
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The  blowing  parameter  can  then  be  expressed  as  follows: 

(a1  + H)  - V a?  - 4a2  + 2a]  H + H2 


1 


B 


Lam 


2a, 


(105a) 


Turb  " a- 


(105b) 


where 

al  - 

0.68  W0-4 

(105c) 

a2  = 

0.08  W0-4 

(105d) 

a3  = 

wo.s 

(105e) 

H = 

AHc/(hr-hw) 

(105f ) 

The  coolant  distribution 

is  then 

given  by 

I 

I 

m = 

B 

(106) 

Tvy  % 

where  qQ  is  the  no  blowing  heat  transfer  to  the  surface  of  the  nose  tip. 
The  cooling  extended  only  to  the  end  of  the  nose  region.  On  the  aft 
cone,  no  cooling  calculations  were  performed. 


4.  TRANSITION  FROM  LAMINAR  TO  TURBULENT 

Any  number  of  methods  for  estimating  the  transition  point  appear 
in  the  literature.  For  simplicity  sake,  transition  was  assumed  to  occur 
when  the  Reynolds  number,  based  on  the  blowing  momentum  thickness, 
equaled  three  hundred,  i.e.. 
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This  value  Is  In  line  with  that  given  by  Van  Driest  (Reference  26)  as 
well  as  Crescl,  MacKenzie,  and  Libby  (Reference  27). 

The  Reynolds  number  based  on  the  blowing  momentum  thickness  if 
given  by 


o u 6** 
^e  e 


(108) 


The  value  for  the  laminar  momentum  thickness  with  blowing  was 
determined  from  Timmer,  Arne,  et  al . (Reference  24)  as 


(109) 


where 


5 **  = 0.4696  S V§-£^- 
o v 3 Re„ 


TRANSVERSE  CURVATURE  AND  VISCOUS  INTERACTION 


(HO) 


The  effect  of  transverse  curvature  Is  to  decrease  the  boundary- 
layer  thickness  and  correspondingly  Increase  the  skin  friction  and  heat 
transfer.  This  effect  Is  small  for  turbulent  boundary  layers  and  hence 
Is  omitted  In  these  calculations.  The  effect  of  boundary  layer  thickness 
on  the  Invlscld  pressure  distribution  is  also  neglected. 
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SECTION  V 


TRAJECTORY  CALCULATIONS 


1.  GOVERNING  EQUATIONS 

The  reentry  trajectory  of  a non-lifting  vehicle  can  be  calculated 
for  a non- rotating  earth  from  the  following  equations: 


1/2  p»  v»  CD  ab  = -m  v»  + m9 sin  y 

m (4>  + y)  ■ mg  cos  y 


(Ilia) 

(111b) 


’G  mV  (<t>  + a) 


172  p«oV®CDAB 


Sketch  14.  Trajectory 

The  centrifugal  force  term,  m $,  can  be  expressed  as 

v®  cos  y ( 

where  R Is  the  radius  of  the  earth  (2.09  x 10+7ft),  and  hr  is  the 
er  « 

altitude  of  the  vehicle.  The  reentry  equations  may  be  written  as 
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! 


then, 

dV  C A 

1ST  = 9 sin  y - 1/2  pw  -D-m£  °°  - (113c) 

and 

&■[?:'  f«er-‘  h6>  ] cos  ■>  <"3,i) 

The  independent  variable  is  the  time,  t,  and  the  unknowns  are  V^,  y.  hg, 
and  X.  The  Runge-Kutta  method  Is  used  to  solve  for  the  above  unknowns 
for  a given  time  step. 


A time  step  of  two  seconds  was  found  to  be  satisfactory  for  altitudes 
above  60,000  feet.  Below  this  altitude  a one  second  time  step  was  used 
and  when  Impact  was  imminent,  the  time  step  was  calculated  from 

At  = 0.9  hg/V^  sin  y (114) 


Impact  was  considered  to  occur  when  the  altitude  was  less  than  10  feet. 

The  Runge-Kutta  method  requires  initial  starting  conditions  which  for  this 
investigation  were  assumed  to  be: 


Altitude 

Velocity 

Angle 

Weight 


300,000  ft 
25,000  ft/sec 
20  deg 
800  lb 


2.  PRESSURE  DRAG 

The  drag  of  the  vehicle  must  be  calculated  at  each  point  in  time 
in  order  to  solve  for  the  unknowns  at  the  next  time  point.  The  total 
drag  of  the  vehicle  is  the  sum  of  the  pressure,  skin  friction  and  base 
drag.  The  skin  friction  and  pressure  drag  are  determined  by  integrating 
the  local  values  over  the  surface  area  of  the  vehicle.  The  pressure  drag 
was  calculated  from 


D 


P 


PJ  sin  9 dA 

00 


(115a) 
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3.  SKIN  FRICTION  DRAG 


The  skin  friction  was  determined  from 


/ T 

•'A 


cos  6 dA 


(115b) 


For  a conical  body  the  element  area  can  be  expressed  as 


dA  = 2iry  ds  * 2iry  hyj  dy 


(116) 


where  the  term  (ds/dy)  for  the  different  nose  shapes  is  given  as  follows: 


Oblate  ellipsoid: 


1 - (y/a)‘ 


(117) 


Flat  nose  - round  shoulder: 

a)  Flat  region:  ds/dy  = 1.0 

ds  rc 

b)  Corner  region:  -r-  = - 

V - <»-v2 

Spherical  arc  - round  shoulder: 


a)  Large  arc: 


(118a) 

(118b) 


(119a) 


ds  c 

b)  Corner  region:  -r-  =*  / — « : 7? 

y y i"c  - (y-y„) 


(119b) 


Aft  cone: 


_J_ 

dy  sin  6, 


(120) 


The  integration  can  be  carried  out  using  any  numerical  integration 
method.  The  Romberg  method,  which  is  explained  in  detail  by  Carnahan, 
et  al.,  (Reference  29)  was  used  for  this  investigation. 
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4.  BASE  DRAG 

The  base  drag  of  a conical  vehicle  was  determined  from  an  empirical 
equation  given  by  Hoerner  (Reference  28). 

c°6  = 5?  02,) 

00 

5.  TOTAL  COOLANT 

Knowing  the  altitude  and  velocity  for  a given  time  the  coolant  was 
calculated  using  the  method  discussed  in  the  preceding  sections.  The 
total  coolant  required  over  the  trajectory  is  of  course  the  summation 
of  the  coolant  used  at  each  time  step.  The  trapezoidal  integration 
method  was  used  to  obtain  this  value  because  the  time  intervals  over 
the  trajectory  were  not  equal  and  the  number  of  points  on  each  trajectory 
varied.  The  expression  used  to  determine  the  total  coolant  was 


Total  coolant  . = Total  coolant  ^ + 1/2  [M^  + ][t. -t^j ] (122) 


AFFDL-TR-78-147 


SECTION  VI 

COMPUTER  PROGRAM  VALIDATION 

The  system  of  equations  presented  In  the  previous  sections  were 
programed  for  use  on  the  CDC  6600  computer,  and  this  section  discusses 
the  validation  of  that  program.  Values  generated  by  the  computer  program 
are  compared  either  with  theory  or  with  data  available  In  the  literature. 

1.  STAGNATION  POINT  VELOCITY  GRADIENT 

A comparison  of  the  stagnation  point  velocity  gradients  with  the 
literature  gives  an  Indication  of  the  validity  of  not  only  the  method  for 
determining  the  velocity  gradient  but  also  the  method  for  calculating  the 
sonic  point  location.  The  stagnation  point  velocity  gradients  for  the 
three  different  nose  shapes  are  plotted  In  Figure  1.  The  open  symbols 
are  from  the  computer  program.  The  closed  symbols  represent  data  from 
the  literature  (References  17  and  34).  As  can  be  seen  the  trends  pre- 
dicted by  the  computer  program  agrees  qualitatively  with  results  In  the 
literature.  The  computer  results  are  somewhat  lower  as  the  flat  face  - 
sharp  corner  shape  Is  approached  because  the  computer  value  Includes  real 
gas  effects  and  because  Vinokur's  hemispherical  value  is  higher  than  the 
value  predicted  by  the  modified  Newtonian  pressure  distribution. 

2.  PRESSURE  DISTRIBUTION 

The  pressure  distributions  over  the  flat  faced  nose  and  the  hemis- 
pherical nose  are  shown  In  Figures  2 and  3,  respectively.  The  pressures 
using  Vinokur's  method,  the  method  developed  In  this  report,  and  the 
modified  Newtonian  method  are  shown  for  comparison.  These  figures  show 
that  the  method  developed  In  this  report  compares  reasonably  well  with 
Vinokur's  results. 

An  overall  indication  of  the  pressure  calculation  method  for  the 
nose  can  be  seen  from  a plot  of  Cn  versus  the  geometric  parameters 
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Figure  1.  Stagnation  Point  Velocity  Gradients 
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b/a,  rc/a,  and  a/RN  given  In  Figure  4.  Also  shown  In  this  figure  are 
the  perfect  gas  range  and  a real  gas  value  for  the  hemisphere.  As  can 
be  seen,  the  computer  program  gives  reasonable  results. 

The  aft  cone  pressure  distribution  for  a hemispherical  nose  tip  is 
plotted  In  Figure  b,  along  with  data  from  the  literature  (References  9 
and  31).  As  expected,  the  values  given  by  the  computer  program  agree 
with  those  from  Krasnov,  (Reference  9)  but  are  somewhat  below  those  of 
Roberts,  Lewis,  and  Reed  (Reference  31).  The  minimum  pressure  value  and 
J location  agree  fairly  well,  but  the  data  of  Roberts,  et  al , is  higher 

beyond  this  point.  The  computer  program  value  for  pressure  on  the  aft 
cone  does  however,  agree  with  the  modified  Newtonian  value. 

3.  NUSSELT  NUMBER  - NO  BLOWING 

The  no-blowing  values  of  the  Nusselt  number  are  plotted  versus  the 
Reynolds  number  in  Figure  6.  The  theoretical  values  are  also  plotted  for 
comparison.  Figure  6 shows  that  the  computer  program  gives  satisfactory 
results  for  both  laminar  and  turbulent  flow.  The  slight  difference 
between  the  computer  values  and  theory  is  due  to  the  diffusion  correction 
in  the  computer  method. 

4.  SKIN  FRICTION  - NO  BLOWING 

The  no-blowing  values  of  the  skin  friction  are  plotted  in  Figure  7 
along  with  the  corresponding  theory.  The  computer  program  gives  slightly 
different  values  because  of  the  diffusion  correction  term  in  the  computer 
program  and  the  assumption  that  C*  = 1 .0  (p*p*/peue)  in  the  theoretical 
calculations . 


The  no-blowing  values  of  the  Reynolds  number  based  on  momentum  thick- 
ness are  shown  in  Figure  8.  This  plot  shows  good  agreement  with  theory 
for  both  the  laminar  and  turbulent  cases. 

5.  COOLANT  DISTRIBUTION 

The  computer  program  is  seen  to  give  satisfactorily  values  for  the 
non  blowing  Nusselt  number,  skin  friction  and  momentum  thickness. 

For  the  blowing  case,  an  overall  indication  of  the  results  of  the  methods 
used  can  be  seen  by  a comparison  of  the  values  from  the  computer  program 
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with  those  of  Schuster,  Gllckman  & Hender  (Reference  32).  For  turbulent 
flow,  Schuster,  et  al . , give  the  following  emperlal  equation  for  the  mass 
flow  distribution. 


5.0  x 10'4  pj^uj-29 
Pr2/3  Rn°-2 


it 


1.0  + 3.03  x 10'8u 


(123) 


where  the  units  are  p(lbm/ft3),  u(ft/sec),  and  RN(ft) 


The  value  for  the  Prandtl  number  Is  0.71,  and  Zj  Is  given  by 

= 8°-6  (0.275)0-8  q2  (124) 

r 0.70 

[1 .0  - 0.168  e^J 

where  0 is  the  angle  with  respect  to  the  axis  of  symmetry  of  the  nose  tip. 


The  mass  flow  is  plotted  versus  the  nondimensional  arc  length  in 
Figure  9.  The  computer  program  gives  reasonably  good  results  with 
slightly  higher  values  for  the  subsonic  region  and  slightly  lower 
values  for  the  supersonic  region. 


6.  TRAJECTORY 


The  validation  of  the  trajectory  section  of  the  computer  program  is 
shown  in  Figure  10,  where  velocity  is  plotted  versus  altitude.  The 
solid  line  is  an  approximation  given  by  Adler  (Reference  33): 


where 


-a  e'ZAh 

V£  e a e 

(125) 

P$L  gSL  CD  S 

2 Zk  W sin  yE 

(126a) 

pj-  In  (p,jE/p);  h >_  75,000  ft 

(126b) 

Zk  - ZM  (75,000-h) ; h < 75,000  ft 

(126c) 

41.5  x 10'6  ft-1 

( 1 26d) 

[13.416  x 10-6  ft'1]2 

(126e) 

60 


I 

i 


Figure  9.  Coolant  Distribution  Over  Hemisphere 
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The  units  of  p are  (lb-sec2/ft4),  g (ft/sec2),  S (ft2),  and  WT  (1b). 
The  equations  assume  that  a and  are  constant,  where  the  computer 
program  allows  for  quasi  variable  CQ  and  y.  Figure  10  shows  that  the 
computer  program  gives  reasonable  results  for  the  trajectory. 
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SECTION  VII 
RESULTS 

The  results  obtained  from  the  Investigation  are  presented  In  this 
section.  For  all  of  the  computer  runs  the  Initial  reentry  conditions 
were  as  follows: 

Altitude  = 300,000  ft 
Velocity  = 25,000  ft/sec 
Trajectory  Angle  * -20  deg 

The  three  different  nose  shapes  were  first  run  with  the  following 
vehicle  conditions: 

Nose  base  radius  (a)  = 1.0  inch 

Length  = 72.0  inch 

Cone  half  angle  = 10  deg 

Weight  = 800  lb 

1 . NOSE  SHAPE  EFFECT  ON  COOLANT 

The  nose  bluntness  of  the  different  noses  was  varied  by  changing 
the  minor  axis  value  (b)  for  the  oblate  nose,  the  shoulder  radius  (r  ) for 

v 

the  flat  face  nose,  and  the  nose  radius  (RN)  for  the  spherical  arc  nose. 
The  results  obtained  are  plotted  in  Figure  11,  which  shows  the  coolant 
required  for  the  trajectory  versus  the  nose  bluntness  parameters,  b/a, 
rc/a,  or  a/R^.  The  spherical  arc  family  is  seen  to  be  the  best  of  the 
families  considered  and  the  flat  face  - sharp  corner  is  the  best  shape, 
requiring  only  about  40%  of  the  amount  of  coolant  of  the  hemisphere.  The 
flat  face  requires  less  coolant  because  of  both  the  smaller  area  and 
lower  overall  heat  loading  of  the  flat  face  as  compared  to  the  hemisphere. 
The  flat  face  surface  area  is  one  half  of  the  surface  area  of  a hemis- 
phere. The  heat  flux  at  the  maximum  heating  conditions  for  the  flat  face 
and  hemisphere  is  plotted  in  Figure  12.  This  figure  shows  that  the  flat 
face  has  a lower  stagnation  point  heating,  that  transition  occurs  later 
for  the  flat  face,  and  that  the  maximum  heating  is  smaller  for  the  flat 
face.  Thus  the  intensity  of  the  heating  for  the  flat  face  is  not  only 
lower  but  also  occurs  over  a smaller  area. 
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As  mentioned.  Figure  11  shows  that  the  spherical  arc  - sharp  corner 
nose  shapes  require  less  coolant  for  bluntness  parameters  less  than  0.70. 
This  is  also  due  to  the  smaller  area  and  heat  loading  of  the  spherical 
arc  nose  shape  as  compared  to  the  flat  face  - round  shoulder  or  oblate 
ellipsoid  shapes.  The  surface  area  for  the  three  different  nose  shapes 
are  plotted  in  Figure  13,  where  the  spherical  arc  shape  is  seen  to  have 
the  least  area  for  a given  axis  ratio.  The  flat  face  - round  shoulder 
has  the  most  area  for  a given  axis  ratio.  This  corresponds  to  the 
coolant  requirements  shown  in  Figure  11.  The  heat  flux  at  the  maximum 
heating  conditions  for  an  axis  ratio  of  0.70  is  plotted  in  Figure  14. 

Even  though  the  spherical  arc  nose  has  a slightly  higher  stagnation 
point  heating,  the  overall  heating  Intensity  is  lower  and  occurs  over  a 
smaller  area. 

For  values  of  a/R^  greater  than  0.70,  Figure  11  shows  the  coolant 
required  for  the  spherical  arc  noses  to  be  the  same  as  the  other  two 
noses.  The  "jump"  at  a/R^  = 0.70  is  due  to  the  "jump"  in  the  sonic  point 
location.  For  a/R^  £ 0.7  (0j  £ 45°),  the  sonic  point  is  assumed  to  be  at 
the  corner.  For  a/RN  > 0.70  ( 0 j < 45°),  the  sonic  roint  location  moves 
toward  the  hemispherical  value.  As  turbulent,  subsonic  flow  has 
associated  with  it  higher  heating  rates,  the  nose  shape  that  has  the 
sonic  point  further  downstream  will  require  less  cooling. 

2.  NOSE  SHAPE  EFFECT  ON  TRAJECTORY 

The  time  to  impact  for  the  data  plotted  in  Figure  11  varied  from 
37.5  to  38.8  seconds.  Practically  speaking  then,  the  time  to  impact  is 
the  same  for  all  of  the  nose  shapes  because  the  nose  drag  was  small 
compared  to  the  total  drag  of  the  vehicle,  7%  for  the  hemisphere  and  11% 
for  the  flat  face.  In  order  to  determine  the  effect  of  the  nose  shape 
on  the  trajectory  as  well  as  the  coolant  required,  the  program  was  rerun 
for  the  oblate  ellipsoid  family  with  larger  nose  base  radii.  The  vehicle 


parameters  were  as  follows: 

a = 2.0"  3.0" 

L = 72.0"  72.0" 

0C  = 9.23°  8.45° 

Rb  = 13.7"  13.7" 

WT  = 1200  lb  1200  lb 
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Figure  14.  Heat  Flux  Distribution  for  Nose  Axis  Ratio  of  0.70 
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For  the  two  inch  nose,  the  nose  contributed  about  37%  of  the  drag  while 
the  three  inch  nose  contributed  about  50%  of  the  total  drag.  The 
coolant  required  is  plotted  versus  (b/a)  in  Figure  15  along  with  the  one 
inch  data.  The  time  to  impact  varied  from  38.0  seconds  to  36.1  seconds 
for  the  two  inch  nose  and  from  41.0  to  36.6  seconds  for  the  three  inch 
nose.  Figure  15  shows  that  the  flat  face  nose  (b/a  = 0)  still  requires 
the  least  amount  of  coolant  for  both  the  two  and  three  inch  nose  radii. 

This  indicates  that  for  hypersonic  impact  the  much  smaller  coolant  flux 
of  the  flat  face  more  than  offsets  the  longer  trajectory  times. 

3.  OPTIMUM  FLAT  FACE  HEIGHT 

The  optimum  nose  shape  determined  from  the  variation  of  parameters 
method  is  in  agreement  with  results  reported  in  the  literature  (References 
1,  2,  and  5).  The  literature  also  shows  that,  for  a given  ratio,  there 
is  a flat  face  height  which  minimizes  the  total  heating  to  the  vehicle. 
These  results  were  obtained  by  assuming  either  a constant  beta  or  constant 
freestream  conditions.  In  reality  the  freestream  conditions  and  beta 
will  change  as  the  flat  face  height  is  changed  for  a given  fineness  ratio. 
To  investigate  the  possibility  of  an  optimum  flat  face  height  without  the 
assumptions  mentioned  above  the  computer  program  was  run  with  the  following 
parameters : 

LT  = 72.0  inch 

Base  radius  = 13.6  inch 

Cone  weight  = 12,000  lb 

Tw  = 2,000  °R 

The  vehicle  weight  was  large  in  order  to  Insure  that  the  high  drag 
vehicles  impacted  hypersonically.  The  flat  face  height  was  varied  so 
that  the  ratio  a/Rg  varied  from  0.1  to  0.9.  The  heat  loading  on  the 
nose,  aft  body  and  the  combined  total  were  calculated  for  the  entire 
trajectory  and  are  plotted  in  Figure  16.  The  beta  of  the  vehicle  at 
maximum  heating  altitude  (Bmh)  is  also  plotted.  Figure  16  shows  that  as 
the  face  height  increases  the  total  heating  to  the  nose  increases  while 
the  aft  heating  and  total  heating  decrease.  The  beta  at  maximum  heating 
altitude  is  seen  to  decrease  significantly.  An  optimum  face  height  is 
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not  observed.  As  the  flat  face  height  increases,  the  drag  of  the  vehicle 
increases  causing  the  velocity  and  hence  heating  intensity  to  decrease. 
The  decrease  in  heating  intensity  over  the  vehicle  more  than  offsets  the 
increase  to  the  nose  because  of  a larger  nose  area.  The  total  coolant 
required  for  the  varying  face  heights  is  shown  in  Figure  17.  As  can  be 
seen,  the  coolant  required  increases  drastically  as  the  flat  face  heights 
is  increased. 
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SECTION  VIII 
CONCLUSIONS 

The  results  of  this  Investigation  can  be  summarized  as  follows. 

For  reentry  vehicles  impacting  hypersonlcally  the  flat  face  - sharp 
corner  nose  is  the  optimum  TCNT  shape.  This  was  found  to  be  true  for 
nose  sizes  that  contributed  up  to  50%  of  the  total  drag  of  the  vehicle. 
Using  a flat  face  nose  shape  reduces  the  amount  of  water  needed,  and 
hence  volume,  about  60%.  This  Is  about  one  pound  savings  in  weight  for 
a one  Inch  nose  radius  and  5.5  pound  saving  for  the  two  inch  nose  radius. 
No  flat  face  height  was  found  that  minimized  the  total  heating  to  the 
entire  reentry  vehicle  for  any  of  the  nose  radii  considered.  If  the 
fineness  ratio  and  weight  of  the  vehicle  are  fixed,  the  overall  heating 
decreases  as  the  flat  face  height  increases  due  to  the  resulting  large 
decrease  in  the  velocity  of  the  vehicle. 

The  effect  of  the  nose  shape  on  the  trajectory  depends  on  the  size 
of  the  nose.  For  low  drag  noses,  7 - 10%  of  total  drag,  the  nose  shape 
has  very  little  effect  on  trajectory  times.  As  the  nose  drag  increases 
to  50%  of  the  total  drag,  the  flat  face  trajectory  time  is  about  12% 
greater  than  that  for  a hemisphere. 

Practically  speaking  it  may  not  be  possible  to  provide  enough 
coolant  at  the  corner  of  the  flat  face.  A rounding  of  the  corner  will 
then  be  necessary.  If  the  corner  can  be  adequately  cooled  for  round 
shoulders  with  rc  less  than  or  equal  to  0.2a,  the  flat  face  - round 
shoulder  shape  is  the  optimum  shape.  If  rc  is  greater  than  0.2a,  then 
the  oblate  ellipsoid  shape  is  the  better  shape. 
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